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Zeeman fields can drive semiconductor quantum wires with strong spin-orbit coupling and in 
proximity to s-wave superconductors into a topological phase which supports end Majorana fermions 
and offers an attractive platform for reahzing topological quantum information processing. Here, 
we investigate how potential disorder affects the topological phase by a combination of analytical 
and numerical approaches. Most prominently, we find that the robustness of the topological phase 
against disorder depends sensitively and non-monotonously on the Zeeman field applied to the wire. 
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I. INTRODUCTION 

Topological quantum information processing promises 
to go a long way towards alleviating the problem of envi- 
ronmental decoherence in quantum computers. In this 
scheme, information is stored and processed by nonlocal 
qubits b ased on quasiparticles with nonabelian quantum 
statistics.'^ These nonlocal qubits are topologically pro- 
tected against local perturbations. Qubit operations are 
based on quasiparticle exchanges which are themselves 
topological in nature and thus insensitive to small varia- 
tions in the operation. 

The simplest type of quasiparticles obeying nonabelian 
quantum statistics are zero-energy Majorana fermions.^ 
Such Majorana bound states are carried by vortices of 
weak-pairing superconductors with spinless p-wave pair- 
ing symmetry, a phase which is widely expected to be 
realized by the v — 5/2 fractional quantum Hall state. ^ 
Starting with the seminal suggestion of Fu and Kane® 
that zero-energy Majorana fermions can in principle be 
engineered to exist in hybrid structures of conventional s- 
wave superconductors and topological insulators, various 
superconductor hybrids have been predicted to support 
Majorana fermions. These hybrids involve, e.g., stan- 
dard two-dimensional electron systems with strong spin- 
orbit coupling,Sl metallic surface states,'^^ semico nduc tor 
quantum wires with strong spin-orbit coupling ,1^^^ or 
half-metallic wires and films .fiSHill 

Realizations based on quantum wires are perhaps par- 
ticularly promising for realizing topological quantum in- 
formation processing since they allow for relatively de- 
tailed scenario^^ of how to manipulate (i.e., create, 
transport, or fuse) Majorana fermions. At the same time, 
disorder is known to often have drastic consequences 
for the electronic properties of one-dimensional electron 
systemd^and for superconducting pairing correlations.'^ 
This motivates us to address the important question 
of how robust these proposed realizations of Majorana 
fermions are against potential disorder. 

The effects of disorder on topological phases which 
support Majorana fermions were studied in relation to 
fermion zero modes inside the vortex core in p-wave su- 



perconductors in the presence of impurities by VolovilPSl 
and, in the one-dimensional case, by Motrunich et al.^ 
Both papers predate the current wave of interest by al- 
most a decade. Motrunich et al. consider a disordered 
p-wave superconductor in one dimension within a renor- 
malization group approach. Their central finding is that 
disorder causes a sharp transition to a non-topological 
phase at a critical disorder strength. More recent work 
on disorder effects on Majorana fermions focuses on the 
strong-disorder limit where the wire breaks up into topo- 
logical and non-topological domains .l^il^Il 

It is the central purpose of the present paper to em- 
phasize that the physics of disorder in proximity-coupled 
semiconductor wires with strong spin-orbit coupling is 
considerably richer. It reduces to the model studied by 
Motrunich et al. (sometimes refered to as Kitaev's toy 
model^P^ only in the limit of a large Zeeman field. In 
contrast, we find that the semiconductor quantum wires 
exhibit additional regimes, depending on the strength of 
the Zeeman coupling, which display characteristic differ- 
ences in the effectiveness of disorder. Most importantly, 
our results can give significant guidance to efforts to re- 
alize Majorana fermions in the laboratory. 

II. SEMICONDUCTOR WIRES 

We consider a semiconductor quantum wire with 
strong Rashba spin-orbit coupling u in an external Zee- 
man field B. For definiteness, we take the spin-orbit 
field and the Zeeman field to point along the x- and z- 
directions, respectively. If the wire is in proximity to an 
s-wave superconductor, the corresponding Bogoliubov-de 
Gennes (BdG) Hamiltonian takes the forniii^ 

'H= (^^+upa^ + Vix)~ f?jT,^ Ba, + AT^. (1) 

Here, p and x denote the momentum and position along 
the wire and A is the proximity-induced pairing poten- 
tial, i.e. the gap in the absence of V, /i, B. The Pauli 
matrices ai (r^) operate in spin (particle-hole) space. Dis- 
order is modeled through a Gaussian white noise poten- 
tial V{x) with zero mean and correlator {V{x)V{x')) = 
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Figure 1. (Color online) Dispersion relations in the absence of disorder for: (a) the quantum wire in the normal phase and 
zero Zeeman field; (b) the quantum wire in the normal phase and finite Zeeman field B <^ Eso, including an illustration (red 
arrows) of a second-order process contributing to the effective disorder potential in the regime of intermediate Zeeman fields; 
(c) the quantum wire in the normal phase and strong Zeeman field B ^ Eso- (d) Quasiparticle excitation spectrum near (full 
line) and at (dashed line) the topological phase transition A = B. The red (dashed) arrows illustrate a second-order process 
to the high-energy subspace near p = 0, which contributes to the random Zeeman field in the low-energy model. 



7(5 (x — x').^'* The disorder strength 7 is related to the 
mean scattering time tq = u/7 at fj, = B — 0. 

Eq. (0 is written in a basis which corresponds to the 
four-component Nambu operator ^' — [i/;^, -i/;^, ^/;|, — i/^l] 
in terms of the electronic field operator 'ip^{x). In this ba- 
sis, the time-reversal operator takes the form T = lUyK, 
where K denotes complex conjugation. The BdG Hamil- 
tonian ([ij obeys the symmetry {H^CT} = 0, with 
C = —iTy. Thus, if is an eigenspinor oiT-L with energy 
E, CT|-0) is an eigenspinor of energy {~E). Zero-energy 
Majorana fermions are characterized by spinors satisfy- 
ing CTIV'o) = \i^o), Majorana modes by CTjV'p) = l^-p)- 



III. REGIMES AND EFFECTIVE 
LOW-ENERGY MODELS 

In the absence of disorder and of the proximity-induced 
gap, i.e., for A — V{x) — 0, the Hamiltonian Eq. ([T|) has 

the particle dispersion E^j^l. = p'^/2m ± [(up)^ + B^]^/^. 
A topological superconducting phase requires that there 
are only two Fermi points, cp. Figs.[l]ja)-(c). In this case, 
the chirality dependence of the spin orientation induced 
by the Ras hba c oupling turns the normal state into a he- 
lical liquid.fii^ The superconducting phase for B < A 
is a conventional, non-topological s-wave superconduc- 
tor. A topological superconducting phase can be real- 
ized when the Zeeman field is larger than A. For strong 
spin-orbit coupling Cgo 3> A (with the spin-orbit energy 
Eso = TOu^/2), it is natural to distinguish three regimes: 



Weak Zeeman fields (ts, 
B - A < A) 



> B > A with 



We can develop a low-energy model by focusing on 
small p where p^ /2m can be neglected relative to the 
spin-orbit term and expanding about the critical point 
S = A at which the p — Q gap closes. The low- 
energy subspace is then spanned by the eigenspinors 
{x\p,±) = (l/2yZ)[l,±l,l,=Fl]'^e^P^ of the BdG Hamil- 
tonian, and the Hamiltonian evaluated in this subspace 
takes the form 



n = 



up — V 
—b + w 



-b + w 
-up + V 



(2) 



Remarkably, all matrix elements of the disorder potential 
in the low-energy subspace vanish so that to this order, 
we find that the disorder fields v — w = 0. Note that 
the eigenspinors of the low-energy subspace correspond 
to Majorana modes. 

The leading contributions to the disorder fields v and 
w are quadratic in the bare disorder potential, emerging 
from virtual excitations of high-energy states shown in 
Fig. [ijd). Adding the contributions from the three rel- 
evant high-energy subspaces near p = and p — ±pf 
{pf = 2mu), we obtain 

v--i:i^ „„"i':^^f,'':\, v.(3) 



v = V 



2A 



^ = "^"^^ TX^"^ + Y ^ , 'To x^V. (4) 



^(p±Pf)" + A2 
A/2 

u'^{p±PFy + A2 



In the low-energy limit, the random gauge field v{x) 
vanishes, while w{x) is a Gaussian white noise Zeeman 
field with {w(x)) = 7/u and {w{x)w{x')) — {w{x))''^ — 
{'-/^ /2uA)S{x - x'). The effective Hamiltonian Eq. ^ 
remains a valid approximation for 7 <^ 7max ^ "A. 



The gap at p = is given by 6 = _B — A in the presence 
of both Zeeman field and proximity-induced supercon- 
ductivity. Thus, as long as |i3 — A| ^ A, the gap at 
p = is much smaller than the proximity-induced gap 
at the Fermi points which is equal to A. In this case, 
the relevant low-energy degrees of freedom for /i ~ are 
those with momenta p near p = 0^ (Note that here we 
take A and B to be positive.) 



B. Intermediate Zeeman fields {eso ^ B ^ A) 

When the Zeeman field is much larger than the 
proximity-induced gap A, but still small compared to 
the spin-orbit energy £50, the gap at p = is of order B 
and thus much larger than the gap at the Fermi points of 
order A. In this case, the relevant low-energy degrees of 
freedom for ~ are those near the Fermi points ±pf- 
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It is important to realize that the spin orientation at the 
Fermi points is dominated by the spin-orbit coupling so 
that electrons at the two Fermi points have almost an- 
tiparallel spins. Thus, the proximity effect is not attenu- 
ated by spin effects, while disorder-induced backscatter- 
ing is strongly suppressed. 

This can be made explicit by projecting the orig- 
inal Hamiltonian Eq. ([ij onto the lower band of 
the normal-state Hamiltonian of the clean wire. 
The corresponding electron and hole eigenspinors are 
\p,e) = [cos(a/2),-sin(a/2),0,0]^e*P^/yt and \p,h) ^ 
[0,0,-sin(a/2),cos(a/2)]^e^P^/A/Z with tana = up/B. 
Evaluating the matrix elements of the full Hamiltonian 
in this low-energy subspace and linearizing the disper- 
sion about the Fermi points, we find a spinless p-wave 
superconductor 

= VcsipX^ ~ Pf)tz + Aoff'Aj^T^ + Vcs{x)XxTz, (5) 

with Fermi momentum pp = 2mu, Fermi velocity Vcs = 
u, and gap function A^ff = A. The disorder potential 
has the correlation function {Vcfi{x)Vcff{x')) — 'jcffS^x — 
x'), corresponding to the scattering time TeS = feff/7cff- 
The Pauli matrices operate in the space of left- and 
right-movers. In Eq. ([s]), we ignore forward scattering by 
disorder, which is not expected to affect the results. 

Unlike for weak Zeeman fields, the effective dis- 
order potential Vcf[{x) has a contribution Vcs{x) = 
{B/4:Sso)V{x) [such that 7cft = {B/4:Sso)'^l\ which is Hn- 
ear in the bare disorder potential. However, this contri- 
bution includes a strong suppression by spin effects and 
for strong enough disorder 7 3> 7iin ^ {Bu){B / EgoY' , 
the dominant contribution to Vcs remains quadratic in 
the bare disorder potential, involving virtual intermedi- 
ate states near p = 0, cp. Fig. [ijb). We find that the 
corresponding magnitude of the effective disorder poten- 
tial is 7eff — j2Bu. Note that for intermediate Zeeman 
fields, the projection to the lower subband remains well 
defined for 7 ^ Bu. 



C. Strong Zeeman fields (B > Eso > A) 

The low-energy theory in this limit has been derived 
previously for the clean case and can be obtained as in 
the case of intermediate Zeeman fields.ESThe relevant de- 
grees of freedom remain those in the vicinity of the Fermi 
points [cp. Fig. [ljc)[, but the spin orientation is now 
dominated by the Zeeman field and electrons at the two 
Fermi points have essentially parallel spin. Consequently, 
the proximity effect is suppressed, while backscattering 
is controlled directly by the bare disorder potential Vix). 
Indeed, the projection [valid for 7 <C Bu(B j 
readily yields a spinless p-wave superconductor Eq. (|5|. 
Unlike in the regimes of weak and intermediate Zee- 
man fields, the detailed parameters differ for the cases 
of fixed chemical potential /i ~ and fixed density 
n — Imuj-K. For fixed /i, one finds pp = {2mBY/'^ . 
v^e = (2B/to)1/^ and Aoff = 2A(e,o/-B)i/2. For fixed 



n, we have pp = 2mu, VcS = u, and A^g = 4:A{eso/ B). 
In both cases, Vcff{x) = V{x) 



IV. PHASE DIAGRAM 

Equation ([s]) is the continuum version of the lattice 
model studied by Motrunich et al.?'^ It has a sharp tran- 
sition to a non-topological phase when AoffTcff = 1/2. At 
the crit ical d isorder strength the density of states v{£) is 
singulailsnillHlZl 



V{e) OC |(WF£/Aeff)ln3(£/Aeff)r\ 



(6) 



whereas for weaker disorder v{e) has a. p ower-law tail 
with an exponent that depends on T^fp^^^, 



v{e) CX e 



4Aotti-ott-3 



(7) 



Although both the topological and the non-topological 
phases have a gapless excitation spectrum in the presence 
of disorder, the topological phase distinguishes itself by 
the presence of Majorana end states at an energy that is 
exponentially small in the wire length L. [221 Alternatively, 
for a wire that is coupled to normal-metal leads, the topo- 
logical phase is characaterized by a negative sign of the 
determinant of the reflection matrix r, whereas det r is 
positive for the nontopological phase.l^Sl 

These results allow us to derive the phase diagram of 
the quantum wire as function of Zeeman coupling and 
disorder in some detail. Consider first the regime dis- 
cussed in Section [III C| of large Zeeman fields at fixed ji, 
for which 1/Teff = 7(m/2B)i/2 and Aoff = 2A{eso/By/^. 
Hence, the phase boundary between the topological and 
non-topological phases is at the critical disorder strength 
7cr = 4Aw. The critical disorder strength decreases with 
the Zeeman field due to the concurring suppression of the 
p-wave gap Aoff when considering the case of fixed den- 
sity, where 1/tcs — j/u and Aoff = 4:A{eso/B) so that 
7cr = 8Au{eso/B). 

In the intermediate regime presented in Section |IIIB[ 
the contribution to Vcg{x) which is linear in V{x) gives 
rise to an elastic scattering rate I/tcs = ^{B /ASsof /u. 
The associated critical disorder strength is 7cr = 
2uA(4eso/S)^. This result is appropriate as long as 
7cr < 7iin, which holds for B > (A/eso)^^^eso- The ex- 
pression for 7cr decreases with the Zeeman field, so that 
the critical disorder strength does not increase with B 
throughout the entire range of B where the phase bound- 
ary between topological and non-topological phases is 
governed by the disorder contribution linear in the bare 
disorder potential. In contrast, the critical disorder 
strength increases with B when B <^ (A/£so)^^^eso, 
where the quadratic contribution to VqS is dominant. 
The quadratic contribution to V^ff gives an elastic scat- 
tering rate of l/roff — ^'^/2Bv?. This yields 7cr — 
2{ABv?Y^'^ for the critical disorder strength. 



In the regime of weak Zeeman fields (c.f. Section III A I, 
the low-energy model in Eq. ^ maps onto a random- 
hopping chain with staggered mean hopping.'^ In this 
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Figure 2. (Color online) Phase diagrams of the quantum wire 
as function of Zeeman field B and disorder strength 7 in a 
log-log plot (a) for fixed density n and (b) for fixed chemical 
potential including numerical results. In (a), the long-dashed 
line delineates the limit of validity of the low-energy models; 
at the short-dashed line, Vcs{x) changes between linear and 
quadratic dependence on the bare disorder V(x). In (b), the 
numerics is based on calculating the determinant of the re- 
fiection matrix for the full Hamiltonian cp. Ref. 1281 7cr,ii 
(dashed line) interpolates between the linear and quadratic 
contributions to 7efi for intermediate Zeeman fields. 



model, an infinite disorder strength is needed to tune to 
the critical poin t if a finite gap is present in the non- 
disordered chain jl^Sinni indicating that any finite disorder 
strength preserves the topological phase. However, we 
still find a dependence of the phase boundary on disor- 
der due to a more basic effect: Since the effective disor- 
der potential is quadratic in the bare disorder V{x), the 
disorder-induced Zeeman field w{x) has not only a ran- 
dom component, but also a finite average (w^x)) = j/u. 
Due to this average, the phase boundary is no longer 
given by 6 = 0, but rather determined by 6 — 7/M = 0, 
or B = A -I- 7/u. The density of states shows a power- 
law dependence on both sides of the topological phase 
transition, with 

i/(e)oce"-\ witha^j-^\4:uA{bu~-/)\. (8) 

Exactly a t the ph ase transition, i'{e) has the Dyson singu- 
larity |6| |25 | 29|30 | ^^th Aeff = j'^/iu^A. Of coursc, these 
conclusions are restricted to disorder strengths smaller 
than the maximal disorder strength 7max ~ uA, where 
the mapping to the Hamiltonian in Eq. ^ is valid. In 
fact, the matching with the regime of intermediate Zee- 
man fields shows that the topological phase is lost once 
disorder becomes of the order of 7max- 

We have also computed the critical disorder strength 
numerically with a method based on a scattering ap- 
proach to the full model ([T]). The method, that is ex- 
plained in detail in Appendix [X] consists in computing 
numerically the reflection matrix r{L) of a wire of length 
L. The existence of the a topological phase (i.e. the pres- 
ence of Majorana bound states) is then signaled by the 
condition det[7-(L)] = —1 [as opposed to det[r(iy)] = 1 
for a non-topological gapped phase] 



V. CONCLUSION 

Our results for the critical disorder strength vs Zee- 
man field are summarized in Fig. 2 showing both the 
analytical and numerical results for the phase bound- 
ary. This phase diagram emphasizes the important im- 
plication that the topological superconducting phase is 
most stable against disorder at the optimal Zeeman field 
Bopt (A/eso)^^^£so- This is valid for fixed chemical po- 
tential [Fig. [2|^b)], but most pronounced for fixed density 
[Fig. [2]Ja)[, where the critical disorder strength 7^ — >■ 
a.s B 00. 

As discussed above, the cause of this effect is the 
vanishing of the effective p-wave gap A^ff in the limit 
B — > cx). It is instructive to compare the large- i? limit 
of the semiconductor model considered here with that of 
a half-metallic ferromagnet, brought into electrical con- 
tact with a superconductor with spin-orbit coupling.'SEIl 
In both models, charge carriers are fully spin polarized. 
Yet, the induced gap Aes is independent of the exchange 
field for the half metals of Refs. [M] and UHl the reason 
being that A^g is set by the strength of the spin-orbit 
coupling in the superconductor. Since Aoff directly de- 
termines the critical disorder strength 7cr, we are thus 
lead to expect that even at fixed density, the critical dis- 
order strength jcr for a semiconductor wire can be made 
to saturate in the limit i? — > 00, if the helical state in the 
semiconductor wire is attributed to spin-orbit coupling 
in the superconductor, rather than spin-orbit coupling 
intrinsic to the semiconductor. 
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Appendix A: Numerical analysis of the topological 
phase 

In this appendix, we describe the algorithm that allows 
us to flnd the boundary between the normal and the topo- 
logical phase in a superconducting quantum wire with 
spin-orbit interaction. The topological phase is marked 
by the existence of bound Majorana end states. Knowl- 
edge of the reflection matrix r(L) of the wire (with length 
L) allows one to determine whether Majorana states ex- 
ist and, thus, to determine the topological properties of 
the wire. Specifically, an odd number of Majorana bound 
state exists precisely if det r — —1.1^ The numerical anal- 
ysis thus focuses on the scattering matrix S{L) of the wire 
and its reflections submatrix r(L). 
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Figure 3. Quantum wire of length L (grey shading) described 
by Hamiltonian H, Eq. (A2|, with unperturbed leads(white) 
described by the Hamiltonian "Ho. The reflection matrix r[L) 
is found by repeatedly concatenating scattering matrices of 
segments with length 5L. 



used previously to study disordered wires with uncon- 
ventional superconductivity^^ and conductivity scaling 
in grapheneP2 In the limit where the segment length 5L 
is small (such that W = Ss — \ deviates only slightly 
from the zero matrix) the scattering matrix of a segment 
can be calculated analytically in linear order in 5L (see 
subsection below). [Note that this approximation is not 
based on the smallness of the perturbations 5% but on 
the smallness of 5L] 

The initial condition S'(O) is to taken to be (open wire) 



The situation we consider is shown in Fig. [3] The scat- 
tering states are chosen as the eigenstates of the Hamilto- 
nian T-Lq of a 'free' wire (lead) with u = 0, A = 0, B = 0, 
but with finite chemical potential /zq, 



"Ho — 



P 

2m 



(Al) 



5 = 



1 
1 



(AlO) 



In order to preserve unitarity to all orders, we set 

Ss~^Ss^{l + ^w\(l- , (All) 



The lead is attached to the central wire segment with 
nonzero it, A, B, and chemical potential fi. The Hamil- 
tonian of the central segment is given by 



where W = Ss — i contains the leading-in-(5L contribu- 
tions. 



(A2) 



1. Scattering matrix of a small segment 



with 

SH = Ar^ + V{x)t^ + upa^T^ - Bcr^ + {fiQ - /i) . 

(A3) 

Here the symbols u, V{x), A, and B are explained in the 
main text. The scattering matrix of a segment of length 
L of the wire has the form 



SiL) 



m r'{L) 
r{L) t'{L) J ' 



(A4) 



To linear order in the segment length 5L the scattering 
matrix Ss can be obtained analytically. The scattering 
states in the leads are indexed by electron-hole charac- 
ter a = ±1, propagation direction (right-moving, left- 
moving) p — ±1, spin in z-direction (up, down) s = ±1, 
and energy measured from the Fermi energy, e. They are 
of the form 



X<7,S 



(A12) 



where t, r and t' , r' are transmission and reflection ma- 
trices for scattering states incoming from the left and the 
from the right, respectively. The scattering matrix of two 
wire segments attached to each other, whose individual 
scattering matrices are 5*1 and ^2, is given by 



'S'l2 = Si® S2 



tl2 ^-12 

ri2 



where 



tl2 

ri2 

''12 



t[[l-r2r[]-%, 
7i+t2[l-r[r2] 



(A5) 



(A6) 
(A7) 
(A8) 
(A9) 



The algorithm to find S{L) consists of repeated con- 
catenations of scattering matrices. The wire with length 
L is split into segments of length SL and the scatter- 
ing matrices Ss = S{SL) of the segments are com- 
bined according to Eq. (A5l. This method has been 



where fc^.f = •\/2m(/io -|- ere) is the momentum, v^,^ — 
k^^/m the velocity, and % is a 4-component spinor in 
Nambu- and spin space. To linear order in SL the contri- 
bution Ss from each term SH can be calculated individ- 
ually. The contributions from the deterministic terms in 
SH are obtained by matching the scattering states 'ipa,p,s,t 
at the interface between the leads and the wire to eigen- 
functions of the wire segment. For the disorder potential 
V{x) the contribution to Ss is random and, to leading 
order in (5L, can be obtained using the first-order Born 
approximation 



^''aa' .ss' 



--(V'a'^-,s',e|V^T,|l^,, + ,,,,) (A13) 



Jo 



-dx 



(A14) 



dx 
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and similarly for t' and r' . From the correlator 
iy{x)V{x')) = ^5{x — x') one finds 



{\5rL',s,?) = (|<5C..f ) ='5,,.,<5.,.,^^7, (A15) 



where vp — Va-,e=o is the Fermi velocity. For each seg- 
ment the contribution of V to Ss is thus taken as a ran- 
dom matrix with zero average and mean square given by 



Eq. (A15 



2. Phase diagram 

Since the wire has a superconducting gap the reflec- 
tion part, r{L), of the scattering matrix becomes unitary 
in the limit of a long wire and for e = the determi- 
nant of r{L) approaches either det ?- ^ +1 or det r — 1 
with the latter indicating the existence of a Majorana 
state. To obtain the phase diagram the algorithm is re- 
peated for different magnetic fields B and different dis- 
order strengths 7 [but with a single disorder realization 
for each disorder strength]. From this the determinant 
det r{L) is calculated and plotted as a function of B and 
7. We have verified that the phase diagram is indepen- 
dent of the choice of the chemical potential fiQ in the 
lead. The data is shown in Fig. 2b together with the 
theoretically predicted phase boundaries. 
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